Numerical study of vortex system quantum melting 
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We report a numerical study of the vortex system in the two dimensional II-type superconduc- 
tors. We have proposed a phenomenological model that takes into account quantum fluctuations of 
Abrikosov's vortices. The results of the quantum Monte-Carlo simulations by the SSE algorithm 
show that the thermal fluctuations are dominated by quantum fluctuations at low temperatures. In 
particular, we demonstrate the possibility of the quantum melting transition for vortex system in 
the temperature region where thermal melting transition is improbable. 
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INTRODUCTION 

Physics of the vortex systems is a very attractive field 
of the modern condensed matter physics due to various 
phases and phase transitions associated with vortex mat- 
ter. Nature of the phase transitions in the vortex system 
is a subject of numerous both theoretical and experimen- 
tal studies (see, for example, the reviews 0,0)- It is well 
known that at low temperatures vortex system is in solid 
phase. Structure (quenched) disorder or other kinds of 
fluctuations have an influence on vortex system state. 
Thus as the thermal fluctuations increase the vortex lat- 
tice melts and system undergoes transition from solid 
to the liquid phase. Some peculiarities of vortex lattice 
melting at the presence of structure defects have been 
obtained by numerical simulations and were discussed in 

Refs. iaiisisEsini- 

Resent experimental results demonstrate an alterna- 
tive mechanism of the order destruction of vortex sys- 
tems. This mechanism is associated by the some re- 
searches with quantum fluctuations in the vortex system 
[3 El El El El So the authors of the paper [H 
investigated magnetic relaxation in MgB 2 thin film and 
made a conclusion on the presence of the quantum creep 
of the vortices at T — > 0. They supposed that the creep 
is induced by quantum fluctuations. The transport prop- 
erties investigations oik- {BEDT - TTF) 2 Cu{NCS) 2 
organic superconductor as well as the measurements 
of dc and ac complex resistivities for thick a — Mo x Si\- x 
films E3 indicated, in the author's opinion, a presence of 
quantum fluctuations too. Note that the quantum fluc- 
tuations are supposed to be essential for thin supercon- 
ducting films and high-layered superconductors, i.e. for 
systems with reduced dimensions. The quantum fluctua- 
tions became potentially significant at low temperatures 
T — > 0, where the thermal fluctuations are negligibly 
small. In other words, quantum fluctuations can give 
rise to melting of the vortex system in the temperature 
range where the vortex system is expected to be in solid 
phase from the classical point of view. 

Though there are many experimental researches of the 



quantum fluctuations influence on properties of the vor- 
tex system, simulations of vortex system with quantum 
fluctuations are practically absent in literature. One ex- 
ception is the paper [l^ where study of two dimensional 
vortex system has been done. However world-line Monte- 
Carlo algorithm in discrete imaginary time used in E3 
did not allow the authors of the paper to simulate the 
processes of creation-annihilation of vortices in external 
magnetic field. 

In this Letter we propose a phenomenological model 
for description of quantum fluctuations in the vortex sys- 
tem and report simulation results which were obtained by 
the loop quantum Monte Carlo algorithm in continuous 
imaginary time. 



MODEL 

To describe the quantum fluctuations of vortex system 
we propose to consider the vortices as hard-core bosons 
and introduce quantum tunneling term into Hamiltonian. 
So, Hamiltonian of vortex system takes the form: 



H=-tJ2 a t a J + h.c. + V ^ 



riiTij 



<hj> 



^hiUi, (1) 



where the first term describes the quantum tunneling 
of the vortices between sites ij with hoping amplitude t 
- the phenomenological parameter for describing of the 
wave function overlap. In the general case vortex can 
hope between any different sites, but for the simplicity of 
the model we take into account the hopping between only 
nearest neighboring sites. The second term describes 
a pairwise interaction between vortices. And the third 
term includes all linear interactions, namely vortex self- 
energy, interaction between vortices and pinning centers, 
interaction between vortices and external magnetic field. 
Quantum behavior of the vortex system is described by 
the first term. The second and third terms are the same 
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FIG. 1: Superfluid density versus temperature for different 
hopping amplitudes t at U = 1.0, h = 0.5. 



FIG. 2: Structure factor Se versus temperature for different 
hopping amplitudes t at U = 1.0, h = 0.5. 



as for classical vortex system (see, for example Ref. |9(). 
In this paper we choose: 



hi = h + fM 



(2) 



(3) 



where h - external field, fii - pinning potential, Ko - 
Bessel function, A - penetration depth, nj - a distance 
between the vortices at sites % and j. 

A study of the system described by equation (JTJ is 
a very complicated problem. Analytic solutions of the 
model ^1 are do not available. Therefore it is neces- 
sary to use simulation methods. The quantum Monte 
Carlo algorithms are the most powerful tools for explor- 
ing quantum many body systems, such as a system l|T[l. 
Quantum Monte Carlo methods allow us to calculate 
thermodynamic characteristics of the model and obtain 
visual qualitative picture of complex physical processes 

ei lis mm 

Note that by use of equation we can describe both 
quantum and classical regimes of the vortex system. 

To analyze the model (QJ we use SSE (stochastic series 
expansion) algorithm modified for long-range potential. 
The algorithm is an exact numerical method and allows 
us to simulate quantum many body system in wide region 
of external fields and temperatur es 1221. The method can 
be implemented in the both site |2lL |22( and interaction 
representations |23j . In the case of long-range potential 



an interaction representation is preferable. For the inter- 
action representation Hamiltonian splits into a diagonal 
part D and a perturbation part V\2?i l24j: 



H = D + V. 
The partition function takes the form: 

Z = E(- 1 )" /„ dri / dr2 

n=0 



(4) 



where 



o JO Jo 
Tr{V{T 1 )V(T 2 )...V(T n )}, 

V{t) = e Tb Ve- rb . 



(5) 



(6) 



After that we choose basis {|a)} where D is diagonal 
and decompose the perturbation part into bond opera- 
tors: 



(7) 



where b - bond index. So we can rewrite the partition 
function in the following way [23j : 



a n=0 T„ 



dT n W(a,T n ,r), 



W = {-l) n {e- pEo e- T ^- E ^){a\ ]J H bp \a),(8) 

p—1 p—1 
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where E p = (a(p)\D\a(p)} , T n a sequence of non- 
diagonal operators. To simplify simulations we add unit 
operators into the operator sequence T n and obtain the 
operator sequence T m . The Monte Carlo simulation is 
carried out by diagonal and loop updates. The simula- 
tion starts with an arbitrary state |a) and the operator 
sequence T m containing only unit operators. The diag- 
onal update contains one attempt to interchange diago- 
nal and unit operators. Note that the diagonal update 
changes the expansion power n by ±1. In a stage of the 
loop update the interchange of diagonal and non-diagonal 
operators is carried out with the hxed expansion power n. 
So the loop update provides the creation-annihilation of 
kinks. At the same time the state of the system \a) can 
be changed. The method allow us to operate directly 
with continuous imaginary time r and grand canonical 
ensemble. 

RESULTS 

We consider a two-dimensional vortex system on trian- 
gular grid with periodic boundary conditions. To analyze 
behavior of the system we calculate a structure factor 

i=l ' 3=1 

and so-called superfluid density 



Superfluid density is introduced in according to the 
conventional determination |25j . Superfluid density de- 
pends on the mean square of winding number that char- 
acterizes a topological configuration of the world lines. 
Note that the nonzero value of p s means the presence of 
the energy dissipation as a result of the vortex quantum 
tunneling, i.e. there is a rise of resistivity. And on the 
contrary as p s — > quantum creep tends to zero and 
the only thermal dynamics of vortices is present (in case 
of current absent). 

Calculations have been done for the vortex system on 
10 x 10 triangular grid with the periodic boundary con- 
ditions. The grid spacing was chosen to be 0.25A. One 
pinning center with the potential \x = 0.1U was intro- 
duced into the system to fix a vortex lattice. Two sets of 
calculations have been done. The temperature is changed 
at the fixed hopping amplitude t and vice verse hopping 
amplitude is changed t at the fixed temperature. 

Fig. ^ shows the superfluid density versus temper- 
ature at different hopping amplitudes t. As seen from 
dependencies at small hopping amplitude t = 0.1 super- 
fluid density is almost absent in whole range of temper- 
atures. While at t = 0.5 well-defined plateau arises on 
graph. The plateau corresponds to the quantum regime 
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FIG. 3: Superfluid density and structure factor S§ versus hop- 
ping amplitude t at temperature T = 0.1C7, U — 1.0, h = 0.5. 

of vortex system. At the same time superfluid density 
drops rapidly as the temperature increase. This corre- 
sponds to transitions of the system to classical regime. 
As mentioned above, the superfluid density has another 
meanning for vortex system then usual one for hard-core 
boson system. Generally if the physical system is in su- 
perfluid state the dissipation is absent. However, as well 
known, high vortex mobility is accompanied by energy 
dissipation and therefore high value of superfluid den- 
sity corresponds to high energy dissipation. So quantum 
fluctuations can give rise to dissipation in the tempera- 
ture range where there is no dissipation in the classical 
approach. 

The structure characteristics of vortex system along 
with the superfluid density are of interest. Figure 2 shows 
the dependencies of the structure factor Sq on tempera- 
ture at different hopping amplitudes t. As seen from the 
curves the structure factor tend to zero at high tempera- 
tures for all considered values of hopping amplitude t. As 
the temperature decreases the structure factor increases. 
Moreover, in the almost classical case t — 0.1 the system 
is more ordered at low temperatures then in the case of 
strong quantum fluctuations t = 0.5. 

To examine a dynamics of the system with the increase 
of the amplitude of quantum fluctuation, we have done 
calculations at fixed temperature T = 0.1U and different 
values of hopping amplitude t. As seen from Fig. [3| an 
enhancement of the hopping amplitude t leads to the in- 
crease of the superfluid density. The last one corresponds 
to the transitions of the system into the quantum regime. 
Simultaneously the structure factor Sq decreases, i.e. the 
vortex system disordering takes place. In other words the 
vortex system melts under quantum fluctuations. 
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SUMMARY 

We have numerically examined the behavior of the 
Abrikosov's vortices system in a two dimensional su- 
perconductor taking into consideration quantum fluctua- 
tions. Simulations have been done in wide range of tem- 
peratures at different intensity of quantum fluctuations. 
We have demonstrated that the enhancement of quan- 
tum fluctuations results in the quantum melting of the 
vortex system. Note that the quantum fluctuations must 
be essential in thin superconducting films, organic su- 
perconductors, high-layered superconductors and can be 
observed in magnetic and transport experiments at low 
temperatures. 
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